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universal results for correlation functions are obtained in the temperature regime T <C A, where A 
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Real time, non-zero temperature (T), correlation func- 
tions of quantum many body systems are directly re- 
lated to many experimentally measurable quantities. For 
strongly interacting systems, there are few quantitative 
results on the relaxation and transport processes that 
are believed to occur at long times at any T ^ 0. Stan- 
dard perturbative methods, as well as many numerical 
approaches, work best in imaginary time, but the ana- 
lytic continuation of such imaginary time results to real 
time is most dangerous, and often fails, in the low fre- 
quency limit. Among the few exact results available are 
those describing relaxational behaviour of an order pa- 
rameter at conformally invariant critical points (obtained 
by exploiting the conformal invariance of the system), or 
by semiclassical methods deep in an ordered state Q, 
both in dimension d — 1. In an important recent de- 
velopment, Altshuler and Tsvelik |2| (AT) obtained real 
time order parameter correlations for the integrable d = 1 
quantum Ising model by the form-factor expansion, and 
their results were in precise agreement with the semiclas- 
sical predictions . 

AT also argued that their results could be extended to 
other gapped, integrable systems, and in particular, to 
the sine-Gordon model; some results for another massive 
integrable model (the 0(3) non- linear sigma model) had 
also been obtained earlier by Fujimoto Q and Konik 
In this paper, we obtain dynamic non-zero T correla- 
tors of gapped one dimensional models in the universal- 
ity class of the d = 1 quantum sine-Gordon field theory 
using a semiclassical method Q : our results are expected 
to apply generic, non-integrable models [f| . Interestingly, 
as we describe in detail below, our semiclassical results 
are qualitatively different from those of AT. One likely 
origin of this difference is that the long time transport 
and spectral Q properties of integrable models, with in- 
finitely many conserved charges [jj , are genuinely distinct 
from those of more generic models; the latter have only 
a small number of conserved charges and are described 
by the semiclassical theory. The generic quantum Ising 
model has only a single conserved charge (the energy den- 
sity) and was possibly too simple a model to expose the 



above distinction. At the end of this paper, we will also 
mention other possible sources of the difference between 
our sine-Gordon results and those of AT and others HQ. 

The imaginary time (r) action of the sine-Gordon 
model is 



A = 
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drdx 
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where $ is the sine-Gordon field, and c has dimensions 
of velocity. It is assumed that the full integrability of the 
model is broken either by additional terms not shown 
above, or by an ultraviolet cutoff that does not preserve 
integrability. Such perturbations are always present in 
any experimentally relevant realization (see Ref. for 
examples of experimental applications), but will be left 
implicit in our discussion below; we will thus freely refer 
to the sine-Gordon model when we mean generic models 
in this universality class. In the integrable sine-Gordon 
theory, g is relevant for 7 < 1, giving rise to a gapped 
phase with particle- like excitations. For l/\/2 < 7 < 1, 
the interactions between these particles are purely repul- 
sive and there are no bound states. It is expected that a 
similar 'purely repulsive' regime exists without full inte- 
grability. 

Our focus here is on the space and time dependent 
correlator of e lJ) * defined in the usual way for arbi- 
trary rj: C${x,t) = ( e *v<£(x,t) e - lv <s>(o,o)^ where $( X)t ) = 

e int ^(x)e^ lUt , H. is the Hamiltonian of the system, and 
the angular brackets indicate average over the equilib- 
rium density matrix of the system. Below, we shall ob- 
tain the space-time dependent C(x,t) for temperatures 
T <C A (A is the gap) in the purely repulsive regime of 
the gapped phase mentioned above. Our answer will be a 
universal scaling function of x, t, A and c, a velocity that 
enters the dispersion relation of the solitonic excitations. 

There are three key observations that allow our exact 
computation for T -C A. The first is that as there is 
an excitation gap, and the lowest lying excitations above 
the gap are 'charged' particles (with charges ±1) that 
represent solitons and anti-solitons in the sine-Gordon 
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field. For T <C A, the particle dispersion e(p) may be set 
equal to its universal low momentum form e(p) = A + 
p 2 c 2 /2A, and they may be treated semiclassicaly using 
Maxwcll-Boltzmann statistics with a velocity probability 
density P(v) ~ e~ Av l 2c ' T . Indeed, their density 
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(1) 



while their r.m.s. velocity v is vt — c{T / A) 1 / 2 . The 
mean inter-particle spacing ~ e +A / T is thus much larger 
than the thermal de-Broglie wavelength ~ c/(AT) _1 / 2 , 
and the particles behave semiclassically .1]. 

The second observation is that collisions between these 
particles are described by their two-particle S-matrix, 
and only a simple universal low-velocity limit of this S'- 
matrix is needed in the low T limit since the particle 
'rapidity' ~ vt/c -C 1. In this limit, in the purely re- 
pulsive regime of the gapped phase, the S'-matrix for the 
process in Fig^]; is generically 
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In other words, the excitations behave like impenetra- 
ble particles which preserve their charge in a collision. 
Energy and momentum conservation in d = 1 require 
that these particles simply exchange momenta across a 
collision (Fig [IJ;). The (—1) factor in Eq. @ can be 
interpreted as the phase-shift of repulsive scattering be- 
tween slowly moving bosons in d — 1. Indeed, the simple 
generic form of Eq. (j2J is due to the slow motion of the 
particles and is not a special feature of relativistic 
continuum theory. It is also worth noting that at special 
free fermion points in the gapped phase (when a map- 
ping to non- interacting fermions exists), the repulsive in- 
teraction between the solitonic particles goes to zero and 
the S- matrix is then perfectly transmitting i. e. it equals 
(— l)6 mi m' fim 2 m' 2 (f° r generic, non-integrable models, ac- 
cessing this free fermion behavior requires tuning an in- 
finite number of couplings, and so this point is not phys- 
ically relevant). The dynamics in this special case is ex- 
pected to be very different from the generic behavior, and 
will be discussed below as a prelude to our results in the 
generic case. 

Finally, we note that in the semiclassical limit at low 
temperatures, for most of the time, the sine-Gordon field 
$ takes one of the values 



$ = 



27m 

7 



(3) 



where n is an integer which increases (decreases) upon 
crossing in the positive x direction a soliton (anti-soliton) 
trajectory. A sample set of soliton trajectories in space- 
time are shown in Fig^ along with the associated values 
of n. Notice the difference in the n values between the 
transmitting soliton free-fermion case, and the generic 
reflecting soliton case described by Eq. 




(c) two-body collision 



~~> x (b) reflecting soiitons 



FIG. 1: A typical set of particle trajectories contributing to 
C<s>(x,t). Each trajectory represents paths moving both for- 
ward and backward in time. The charge moving on a tra- 
jectory is indicated by different colors (color online): Red 
corresponds to soiitons (m = +1) and blue to anti-solitons 
(ra = — 1). The numbers in the domains are the local values 
of the integer n in Eq. 10 



We now present our final results for C$(x, t), deferring 
their derivation till later. For the transmitting soliton 
case in Fig^,, an elementary computation shows that 

C$(x,t) = Aexp (-2,0 sin 2 (7^/7) (\x - vt\)) , (4) 

where the angular brackets represent an average of v over 
P(v) and A is an r\ and 7 dependent amplitude related 
to the vacuum expectation value of e 4 ' 7 *. Remarkably, 
Eq. is precisely the result obtained by AT for the 
integrable sine-Gordon model for generic values of 7 even 
away from the free-fermion point. So it appears that the 
dynamics of the integrable model is rather like that of 
non-interacting particles. 

For generic reflecting soliton case in Fig^, with the 
S-matrix in Eq. J5J (which we claim is always the ex- 
perimentally relevant case), we obtained a very different 
result. We found 



U (2iq r Q, 2i^J q r q t ) 



U Q (2iqiQ, 2iyjq r qi) - iU\(2iq r Q, 2iy/q r qi) 



-iU^iqiQ, 2iyjq r qi) - I (2y/q r qi) 



where 6 = cos(27r?7/7), q r = p J^j^ dvP(v)(x — vt), qi 
p dvP(v)(x — vt), Iq is the modified Bessel function, 
and U0.1 are the Lommel functions of two variables 0. 
(Note that all x and and t dependence is through q r and 
qi which can both be written as functions of the scaling 

variables x — px = x/^ x and i — cpJ ^is^ = 



(5) 
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This remarkable expression is qualitatively different in 
its asymptotic behavior for large t and x from the cor- 
responding result at the free fermion point. One way to 
see this is to obtain the correlator of the soliton charge 
density p = {<y/2n)d x §, C m {x,t) = (p(x,t)p(0,0)), by 
expanding Eq. JSJ to order rf and taking two x deriva- 
tives. The result is identical in form to the diffusive corre- 
lation function obtained for the conserved (spin) density 
m RefH using the analysis of Refs H E3— here, the dif- 
fusion constant has value D m = j^e A ^ T ■ A more direct 
indication of diffusive behavior follows from an asymp- 
totic analysis of our result. For instance, for large t and 
x ~ using q r /i ~ t/£ t ± x/2^ x and standard results 
for the asymptotics of Bessel and Lommel functions 0, 
we find: 



C$>(x, f) 



A 



i + e 



-x 2 /(4D 4 ,t) 



p(l -Q)J (AttD^) 1 / 2 



(6) 



where the diffusion constant £)$ = c\[T j p\p2ixK — 
■§]±e A ^ T is identical to the diffusion constant D m that 
describes the diffusive behaviour of the soliton charge cor- 
relator. This may be interpreted by noting that the field 
values t) and <&(0, 0) are perfectly correlated if the 
space-time points (0,0) and (x,t) lie in the same domain, 
and this happens when the domain walls bounding this 
domain diffuse from the neighbourhood of the origin to 
the space-time point (x, t); the equality of the two diffu- 
sion constants then suggests that such events dominate in 
the average over all space-time histories. In this context, 
it is however important to bear the following subtlety in 
mind: The asymptotic expression Eq. © is only valid 
as long as |0| < 1 — x£ t / (2t£ x ) , and as such cannot be 
used directly to expand in r\ about r] — (which of course 
corresponds to 9 = 1). Put another way, the long time 
limit does not commute with the — > 1 limit. 

We now describe the derivation of Eqs. and JSl. 
We represent C(x,t) as a 'double time' path integral 1], 
with the e~ lHt factor generating trajectories that move 
forward in time for each particle, and the e iWt producing 
trajectories that move backward in time. In the classical 
limit, stationary phase is achieved when the trajectories 
are time-reversed pairs of classical paths (Fig^l. Each 
trajectory has a charge label which obeys Eq. at each 
collision; however as each collision contributes both to 
the forward and backward trajectories, the net numeri- 
cal factor is simply +1. All of this implies that the 
straight lines (rays) y M (i) = y M + v^t (p = 1,2, ...M 
)in Fig ^ are independently distributed with a uniform 
distribution in space for the intercepts y^, and with in- 
verse slope determined by the velocities which are 
independently distributed according to the probability 
density P(v). The charge, m, is assigned randomly to 
each trajectory at some initial time t = with probabil- 
ity Pm/p — 1/2, but then evolves in time as discussed 
above (Fig^l. We have assumed a large system size L, 



and will eventually take the limit L — > oo, M — > oo with 
total density p = M/L fixed and given by Eq. Q . Note 
that the charge travelling on a ray changes in time when 
the iS-matrix is reflecting, and charges on rays are mu- 
tually uncorrelated only at t = — it is the charges rrik 
(A; = 1,2,... M labelled starting from the left) travelling 
on the complicated zig-zag trajectories Xk (t) that remain 
uncorrelated and constant in time. Conversely, when the 
S'-matrix is perfectly transmitting at special free-fermion 
points in the phase diagram, charges on rays remain un- 
correlated and constant in time (since the rays and the 
trajectories coincide in this case). In both cases, the field 
<fr(x,t) is given in terms of soliton trajectories by 



$(x, t) = — m k 0(x - x k (t)). 
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(7) 
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To proceed further in the generic reflecting case, we 
note that the spatial sequence of domains encountered 
at any given time as we move from x = — oo to x = +oo 
is invariant under the time evolution. Therefore, it is 
valid to label the sequence of domains from the left by 
an index I whereby the field $ takes value <!>; indepen- 
dent of t for all x in the I th domain. In doing so, we 
resolve the singularity associated with each collision by 
stipulating that a zero length domain of the appropriate 
field value continues to exist at the instant of the colli- 
sion. Now, let space-time points (x, t) and (0, 0) lie in 
the 1% and l\ h domain respectively — such an assignment 
is well-defined precisely because the sequence of domains 
is time- invariant. The difference Q(x,t) — $(0,0) is thus 
equal to $; 2 — $/ x . To obtain the latter, it is extremely 
useful to work at t = and write 



4>, 



4i 



2?r 

1 



i=i 



(8) 



where mi,,,, are the charges of the s = \h — h\ domain- 
walls (solitons or antisolitions) encountered when going 
from domain I2 to domain l\ at t — 0. The advantage 
of this procedure is clear: The ?7ii...|i 2 _i 1 | are uncorre- 
lated independently random variables that take values 
±1 with equal probability, and as such, may be averaged 
over independently of one another for fixed s. 

The number s is of course a complicated function of 
x, t and all the initial positions y^ and velocities of 
all the rays. However, since the t = soliton charge 
on a ray is not correlated with the y^ and v^, we may 
perform the average over these s soliton charges first for 
a given realization of the y M and , and then average the 
resulting expression over all the y^ and v^. The result 



,ii7*(x,t)„-ii7*(0,0)\ 



/ charge 



of this charge averaging is ( 

Ae s 

Next, we note that s may be obtained by moving from 
the spacetime point (0, 0) to point (x, t) along the straight 
line joining the two, and keeping track of the numbers 
k r and ki of rays that cross this straight line from the 
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right and left respectively. After every such right inter- 
section, we move one domain to the right, while every 
left intersection results in moving back one domain to 
the left — thus, s is precisely equal to \k r — Since 
initial positions and velocities of the rays are indepen- 
dent random variables, the probability of having k r right 
intersections and k\ left intersections may be calculated 
straightforwardly by elementary means in terms of q r , 
the probability that the straight line connecting point 
(0, 0) to point (a;, i) has a right intersection with a given 
ray, and qi, the corresponding probability for left inter- 
section. This, in conjunction with the result for charge 
averaging above, then allows us to write the following 
combinatorial expression 



rescaling of integration variables then gives: 



C$(x,i) = A 



lim 

M^oo 



M-k r 



k r , ki 



M\q^qf'(l-q r -q,) 
k r \ki\{M - k r - h)\ 



. |fer-fc;| j (9) 



where the primed sum is over all k r and ki ranging from 
to M with the constraint that k r + ki < M, q^ r — qi, r /M, 
and we have assumed without loss of generality that x 
and t are positive. 

Before we proceed further, we note that the perfectly 
transmitting case is considerably simpler: Since the tra- 
jectories coincide with rays, and charges on rays are in- 
variant in time and uncorrelated, &(x,t) — <&(0, 0) can 
be obtained by moving along the dotted line connecting 
(0, 0) and (x, i) and counting the total number k r + ki of 
rays intersecting this line (this is the number of uncor- 
related random charges we cross in going from (0, 0) to 
{x, t)) and multiplying by 2ir/~f. The rest of the argument 
above is unchanged and we obtain a simpler combinato- 
rial expression with the \k r — fc;| in the exponent of the 
cosine in Eq. 10 replaced by k r + fcj. Taking the ther- 
modynamic limit using the trinomial formula gives the 
result Eq. (gj. 

Returning to the reflecting case, further progress re- 
lies crucially on a device suggested to us by D. Dhar. 
The artifice involved consists of viewing the sum of a 
subset of terms with fixed k r — fc; = n as the coeffi- 
cient of z n in the Laurent expansion of some function 
of a complex variable z. For n > 0, this is achieved by 
identifying the corresponding partial sum with the co- 
efficient of z n in the trinomial expansion of Jm{®z) = 
(1 — q r — qi + Qzq r + qi/Qz) K1 . Similarly, for n < 0, the 
identification is with the coefficient of z n in the trinomial 
expansion of /m(z/Q). The advantage of this approach 
is now clear: Each Laurent coefficient can be written as 
a contour integral around a contour enclosing the origin, 
and the two sums over n for either sign of n are simply 
geometric series which can be evaluated in closed form. 
Moreover, the thermodynamic limit can be taken explic- 
itly for both sums upon choosing in each case a contour 
along which the geometric series converges. A convenient 



C$(x,<) = A 




q r qi{w+w ) 



(10) 



Here C a is an anti-clockwise contour enclosing the origin 
with \w\ > \Q\y/q r /qi and Cj, goes clockwise around the 
origin with \w\ < \J q r /qi/\<d\. We now expand the ex- 
ponentials in a Laurent series and integrate each term 
to finally arrive at Eq. JSJ upon resumming the resulting 
series. 

The main issue raised by our results is that of the ori- 
gin of the distinction between the semiclassical and form- 
factor results. Both methods make subtle unproven as- 
sumptions about orders of limits: the commuting of the 
t — > oo limit with either the semiclassical (in which the 
low momentum limit is taken in obtaining Eq. J5J) or the 
form factor expansions. It could be that such an assump- 
tion is invalid in one of the methods, and only the other 
result is generically valid in the long-time limit. However, 
we suspect that the distinction reflects a fundamental 
physical difference between generic and integrable sys- 
tems, and the two results apply in their respective cases. 

After this work was completed, we became aware of 
independent work by Rapp and Zarand[ll] who find that 
diffusive behavior appears in the correlators of the one- 
dimensional Q-state quantum potts model for Q > 2. 
Remarkably, the correlation functions in the two cases 
differ only in the choice of the parameter which takes 
the value <d = —1/{Q— 1) in the Potts case. 
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